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o : 

D . Abstract. We present the first models of Saturn and Jupiter to couple their evolution to both 

' a radiative-atmosphere grid and to high-pressure phase diagrams of hydrogen with helium. The 

purpose of these models is to quantify the evolutionary effects of helium phase separation in Saturn's 
deep interior. We find that prior calculated phase diagrams in which Saturn's interior reaches a 
region of predicted helium immiscibility do not allow enough energy release to prolong Saturn's 
, cooling to its known age and effective temperature. We explore modifications to published phase 

^ ■ diagrams that would lead to greater energy release, and find a modified H-He phase diagram that is 

^sO \ physically reasonable, leads to the correct extension of Saturn's cooling, and predicts an atmospheric 

OO . helium mass fraction Yatmos in agreement with recent estimates. We then expand our inhomogeneous 

CO ' evolutionary models to show that hypothetical extrasolar giant planets in the 0. 15 to 3.0 Jupiter mass 

range may have TeffS, 10-15 K greater than one would predict with models that do not incorporate 
. helium phase separation. 

o \ 

'p,: INTRODUCTION 

o ■ 

^ \ The interiors of Jupiter and Saturn, extrasolar giant planets (EGPs), and brown dwarfs 

^ ■ (BDs) are all described by similar physics: these bodies are mainly composed of liquid 

. . ! metallic hydrogen, and their interior energy transport mainly occurs through efficient 

\ convection [T], leading to largely isentropic interiors. Jupiter and Saturn, whose radius, 

^ ■ mass, luminosity, and age are known precisely, can serve as calibrators of thermal- 

. history calculations for the entire class of objects. They can provide a test of the ad- 

equacy of the diverse physical models, including interior thermodynamics, heat trans- 
port mechanisms, and model-atmosphere grid, that enter into the general thermal-history 
theory for EGPs and BDs. However, at very low effective temperatures (~ 100 K), the 
corresponding interior temperatures may become low enough for phase separation of 
abundant interior components to occur, and this effect must be quantitatively evaluated 
before Jupiter and Saturn can be used as calibrators. This work provides a quantitative 
assessment of inhomogeneous evolution in Jupiter, Saturn, and low-mass EGPs. 
Relevant calculations and data for Jupiter and Saturn follow: 

• Saturn is currently over 50% more luminous than one would predict with a homo- 
geneous, adiabatic cooling model. Saturn models reach the planet's known T^ff of 
95.0 K in only 2.0 to 2.5 Gyr Jl. 

• For Jupiter, homogeneous, adiabatic cooling models allow Jupiter to reach its 
known Tgff of 124.4 K in ~ 4.7 Gyr. This is a good match. 
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FIGURE 1. Temperature-pressure plot of the interiors of Jupiter and Saturn at f = 4.56 Gyr, superim- 
posed on a hydrogen phase diagram (see text for details). The upper boundary of the horizontally-hatched 
region shows the minimum temperature at which He is fully miscible in metallic hydrogen with a mass 
fraction Y = 0.27, while the lower boundary shows the minimum temperature corresponding toY = 0.21, 
according to HDW 0| theory. The lower vertically-hatched region shows the same He miscibility limits 
according to Pfaffenzeller, et al. ff\, while the upper vertically-hatched region shows the modification to 
the Pfaffenzeller et al. theory that gives a realistic prolongation of Saturn's age. With this phase diagram, 
when He rains down, it falls all the way down to the core. Liquid metallic hydrogen (H+), liquid molecular 
hydrogen (H2), and the maximum pressures and temperatures reached by laser shock data are also shown. 



• The atmospheres of both Jupiter and Saturn are depleted in helium relative to the 
expected protosolar helium mass fraction of ~ 0.27. Yatmos for Jupiter is 0.231 10] 
and for Saturn is likely 0.18-0.25 |0], but is poorly constrained. 

HYDROGEN AND HELIUM UNDER HIGH PRESSURE 

Under extreme pressure dense molecular hydrogen dissociates and ionizes to become 
liquid metallic hydrogen. The transition is likely continuous over a pressure of 1-5 Mbar 
at the temperatures of interest (~ 10^ K) in giant planets. Also at these temperatures, 
helium, which makes up about 8% of the atoms in a solar composition mixture, likely 
becomes immiscible in liquid metallic hydrogen (3,^,^. Figure [l] shows the interior 
temperature-pressure profile of Jupiter and Saturn superimposed on a hydrogen phase 
diagram. 

Detailed calculations on the dynamics and distribution of helium in giant planets have 
been performed [ 8] . These authors found that when helium becomes immiscible in liquid 
metallic hydrogen, the composition that separates out is essentially pure helium, and this 
helium on fairly short timescales (relative to the convective timescale) will coalesce to 
form helium droplets. These droplets, being denser than the surrounding liquid metallic 
hydrogen, will fall through the planet's gravitational field. If the droplets reach a region 
where helium is again miscible at higher concentration, they will redissolve, enriching 
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FIGURE 2. Evolutionary models including the phase separation of helium from liquid metallic hy- 
drogen. The dotted curve for Saturn and the solid curve for Jupiter are for homogeneous evolution. The 
dashed curve includes the phase diagram of HDW|6|, while the solid curve for Saturn uses the proposed 
modified phase diagram (MPfaff). The modified phase diagram allows both Jupiter and Saturn to reach 
their known ages and Tef/s. The interior of Jupiter begins He phase separation at f « 5 Gyr. 



the deeper regions of the planet in helium. Helium would be lost from all regions 
with pressures lower than the pressures in the immiscibility region, since the planet 
is fully convective up to the atmosphere. Excess helium would be mixed down to 
the immiscibility region and be lost to deeper layers. This "helium rain" could be a 
substantial additional energy source for giant planets. 



CALCULATIONS 

We find that the calculated phase diagram of HDW [6], which is essentially equivalent 
to that of Stevenson [5], are inapplicable to the interiors of giant planets, if helium 
phase separation is Saturn's only additional energy source. As Figure |2lshows, this phase 
diagram prolongs Saturn's cooling 0.8 Gyr, even in the most favorable circumstance that 
all energy liberated is available to be radiated, and does not instead go into heating the 
planet's deep interior. As we show in our published work [2], we find that a modified 
version of the phase diagram of Pfaffenzeller et al. with a higher temperature for the 
onset of helium immiscibility (see Figure [T]), allows Saturn to reach its known Tgff and 
age, while Jupiter evolves homogeneously until t ^ 5 Gyr. Saturn's Yatmos falls to 0.185 
at 4.5 Gyr, which is at the low end of Saturn's derived value [4]. This modified phase 
diagram can be applied to various hypothetical giant planets |9]. We follow the evolution 
of planets with masses from 0.15 Mj (1/2 Saturn's mass) to 3.0 Mj. Figure |3] shows the 
evolution of these planets with and without the effects of helium phase separation. These 
planets are in isolation and possess 10 Mgarth cores. At Gyr ages the model planets 
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FIGURE 3. Evolution of the Teff for planets of mass 0.15 to 3.0 Mj for our standard models with no 

stellar irradiation and 10 Mgarth cores. The dotted lines are models without helium phase separation, while 
the solid lines include the effects of helium phase separation on the planets' cooling. The top curve is for 
the highest mass planet while the bottom curve is for the lowest mass planet. 



undergoing phase separation can have TgffS 10-15 K higher than the homogeneous 
models, making the planets somewhat easier to detect. 
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